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ABSTRACT 


A  computer  program  for  calculating  plane  wave  transmission  through 
an  inhomogeneous  cold  plasma  slab  is  described.  The  program  accepts 
arbitrary  profiles  (including  jumps)  in  number  density  and  collision  frequency. 
It  prints  out  the  electric  and  magnetic  field  profiles  for  the  angles  of  incidence, 
direction  of  incidence,  and  polarisations  of  the  wave  that  are  called  for.  It 
also  gives  the  reflection  and  transmission  coefficients  in  each  case. 
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I.  INTRODUCTION 


In  estimating  the  attenuation  and  reflection  of  plane  electromagnetic 
waves  by  a  plasma  layer  it  is  often  adequate  to  replace  the  layer  by  an 
equivalent  one  with  uniform  properties,  for  which  the  transmission  loss  is 
given  by  a  simple  formula  (Ref.  1).  On  the  other  hand,  in  cases  where  the 
plasma  properties  vary  significantly  in  the  space  of  a  wavelength,  the  uniform 
slab  approximation  is  not  justified.  Furthermore,  in  studying  breakdown 
under  intense  electromagnetic  fields  one  is  interested  in  the  electric  field 
profile  in  the  vicinity  of  points  where  the  dielectric  constant  goes  to  zero, 
and  then  it  is  clearly  necessary  to  take  account  of  the  variation  in  the 
dielectric  constant.  Finally  it  is  useful  to  have  a  systematic  way  for  com¬ 
paring  existing  programs  (Refs.  2  through  4)  and  for  checking  estimates 
based  on  equivalent  homogeneous  slabs  against  exact  calculations. 

For  these  reasons  it  was  decided  to  write  a  computer  program  for 
the  inhomogeneous  case  that  would  be  as  flexible  and  easy  to  use  as  possible. 
The  present  report  describes  the  program  and  gives  a  listing  of  it  in  the 
hope  that  it  may  also  be  of  value  to  others. 

The  program  applies  to  a  cold  plasma  that  is  characterized  by  its 
electron  number  density  and  electron  -neutral  collision  frequency  profiles. 

It  is  designed  to  be  as  flexible  as  possible  in  the  following  respects.  The 
user  is  free  to  input  any  arbitrary  profiles  of  number  density  and  collision 
frequency  (including  discontinuities)  by  listing  the  values  of  these  quantities 
over  any  arbitrary  set  of  corresponding  abscissa  values.  He  can  call  for 
the  calculations  to  be  performed  for  any  number  of  frequencies  and  for  any 
uniformly  spaced  sequence  of  angles  of  incidence,  and  for  either  or  both 
polarizations  of  the  incident  wave.  Finally  he  can  have  the  wave  be  incident 
on  the  slab  from  either  side.  The  computer  output  lists  the  collision  frequency 
and  number  density  profiles  that  were  input.  It  then  tabulates  the  transverse 
electric  and  magnetic  field  profiles,  the  dielectric  constant  profile,  and  the 
profile  of  the  absolute  value  of  the  total  electric  field  through  the  slab  for 
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each  case  (angle  of  incidence,  frequency,  and  polarization)  that  was  called 
for.  Finally  it  gives  the  reflection  and  transmission  coefficients. 

The  next  section  shows  some  of  the  details  of  the  formulation.  The 
following  one  gives  the  exact  format  in  which  the  data  must  be  input  and 
shows  a  sample  of  the  output  in  an  illustrative  case.  Finally  the  program 
itself  is  listed  in  the  Appendix. 
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FORMULATION 


The  plasma  slab  is  considered  to  be  contained  in  the  region 
0  <  X  <  d  with  the  xy  plane  being  the  plane  of  incidence.  For  a  plane 
wave  incident  on  the  slab  at  an  amgle  6  with  respect  to  the  x-axis,  Max¬ 
well's  equations  take  the  form 


dE 
_ z 

dx 


=  -iuB 

y 


-ioj 


k2(p2.K)E^ 


=  (p/c)E^ 


for  the  TE  mode  (for  which  E^, 


and  B_  are  zero)  and 
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for  the  TM  mode  (for  which  B  ,  B  ,  and  E  are  zero).  Here  B  =  sin6 

X  y  z  * 

and  kg  =  w/c,  where  w/2w  is  the  frequency.  The  dielectric  constant  K  for 
a  cold  plasma  is  given  by 

2 

Wp 
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where  is  the  electron-neutral  collision  frequency  and  where  b>^  , 
the  square  of  the  plasma  frequency,  is  proportional  to  the  electron 
number  density  n.  If  n  is  measured  in  electrons /cm  then 

=  3.  186  X  10^ n 
P 

For  the  TE  mode  the  computer  gives  the  profiles  of  E^,  B„,  and 
|e^|  .  For  the  TM  mode  it  gives  the  profiles  of  E^,  and  V|Ej^+  (e|^  . 
The  electric  field  components  E^  or  E^  are  always  normalized  to  unity  on 
the  incident  side.  The  K  profile  is  given  in  both  cases. 

Outside  the  plasma  slab  the  waves  are  taken  to  be  of  the  form 

(1-R)  E  =  exp(ik0CO8dx)  -  R  exp(-ikQCOsex)  for  x  <  0 

(1-R)  E  =  T  exp/ik^cosS  x)  for  x  >  d 

where  E  is  any  component  of  the  electric  field.  Solving  for  R  and  T  we 
find 

E^'  -  ikQCose  E| 

E'  -f  ikgCOsS  Ej 

(l-RXE^/E^)  expi-ik^cose  d) 

in  terms  of  Ej  and  E2.  which  are  the  values  of  E  at  x  =  0  and  x  =  d, 
respectively.  (The  prime  denotes  the  x  derivative.)  The  computer  furnishes 
these  values  together  with  their  decibel  equivaiants  (defined  as 
20  logjQlRj  and  20  iogj^jTl ,  respectively). 


R  = 


T  = 
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m.  DATA  INPUT 


The  data  are  input  on  regular  80-coluinn  Fortran  cards.  At  least 
7  cards  are  required  for  each  case.  A  particular  case  consists  of  a 
plasma  slab  defined  by  its  thickness  and  its  profiles  of  number  density 
and  collision  frequency  together  with  a  set  of  plane  waves  having  1  to  10 
different  frequencies  and  any  number  of  angles  of  incidence.  Any  number 
of  cases  can  be  run  by  stacking  the  corresponding  groups  of  cards.  The 
cards  required  to  represent  one  case  are  as  follows: 


CARD  1  contains  6  two-digit  numbers  placed  in  first  12  columns  as  follows; 
(A  one-digit  number  is  preceded  by  0,  e.  g.  ,  3  becomes  03.  ) 


Columns  1,2:  Number  of  stations  at  which  number  density  is  given. 

(See  CARD  2.  )  Maximum  allowed  is  30. 

Columns  3,4:  Number  of  stations  at  which  collision  frequency  is  given. 

(See  CARD  3.  )  Maximum  adlowed  is  30. 


Columns  5,  6: 
Columns  7,8: 
Columns  9.  10: 

Columns  11,  12 : 


Number  of  frequencies.  Maximum  allowed  is  10. 

Number  of  angles  of  incidence.  (See  CARD  6.  ) 

Mode;  01  for  TE,  02  for  TM  and  00  (or  blank)  for  both. 
(Only  one  mode  is  computed  at  normal  incidence.  ) 

Direction  of  x-component  of  incident  wave:  00  (or  blank) 
if  to  the  right  (i.  e.  ,  the  direction  in  which  the  data  on 
CARDS  2  -  5  is  entered),  01  if  to  the  left. 


CARD  2  lists  stations  (values  of  x/d)  at  which  number  densities  are  to  be  given. 
These  values  must  be  listed  in  increasing  order  starting  with  0. (representing 
the  front  of  the  slab)  and  ending  with  1.  (representing  the  back).  Each  one 
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(including  the  last)  requires  a  decimal  point  and  occupies  5  columns,  so 
the  total  number  of  columns  is  5  times  the  first  number  entered  on  CARD  1. 
If  more  than  80  columns  are  required,  two  cards  must  be  used,  and  the 
second  then  counts  as  the  continuation  of  CARD  2.  If  at  any  station  two 
values  of  the  number  density  are  to  be  given  (corresponding  to  a  discontinuity 
in  the  profile)  then  the  corresponding  station  must  appear  twice.  Thus  a 
homogeneous  slab  requires  4  stations  viz.  ,  0.  ,  0.  ,  1.  ,  1.  . 

CARD  3  lists  stations  at  which  collision  frequencies  are  to  be  given. 
Remarks  under  CARD  2  apply. 

9  3 

CARD  4  lists  number  densities  in  multiples  of  10  /cm  corresponding  to 
the  stations  listed  on  CARD  2.  The  first  and  last  values  correspond  to  the 
media  carrying  the  incident  and  transmitted  waves  respectively,  and  will 
normally  be  0  (for  vacuum).  Each  value  occupies  5  columns,  so  the  total 
number  of  columns  used  is  the  same  as  on  CARD  2.  The  smallest  possible 
number  density  is  10^/cm^  (entered  as  .  0001)  and  the  largest  is  10^^/cm^ 
(entered  as  99999).  Values  between  10‘^/cm'*  and  lO^’/cm'’  are  entered 
as  5  digit  numbers  without  a  decimal  point. 

9-1 

CARD  5  lists  collision  frequencies  in  multiples  of  10  sec  corresponding 
to  the  stations  listed  on  CARD  3.  Remarks  under  CARD  4  apply. 

CARD  6  contains  3  numbers,  each  occupying  5  columns,  including  the 
decimal  point. 

Columns  1-5:  Slab  thickness  in  centimeters. 

Columns  6-10:  First  angle  of  incidence  to  be  run;  0  (or  blank) 

corresponds  to  normal  incidence.  Thereafter  this 
value  is  incremented  by  the  number  appearing  in 
Columns  11-15  of  this  card  until  a  total  number 
of  angles  of  incidence  has  been  run  equal  to  the 
number  appearing  in  Columns  7,8  of  CARD  1. 
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Columns  11-15: 


Interval  in  degrees  between  successive  angles 
of  incidence. 


CARD  7  lists  wave  frequencies  in  GHz.  Each  value  occupies  5  columns, 
including  the  decimal  point.  The  total  number  of  frequencies  listed  must 
not  be  less  than  the  number  appearing  in  Columns  5  and  6  of  CARD  1. 

As  a  particular  (rather  artificial)  example  we  consider  a  plasma 

slab  1.  50  cm  thick  and  having  an  electron  number  density  which  starts  at 

10^^/cm^  at  x/d  =  0,  falls  linearly  to  lO^Vcm^  at  x/d  =  .63,  at  which 

point  it  jumps  to  10  /cm  ,  and  then  fails  off  linearly  to  zero  at  x/d  =  1. 

We  take  the  collision  frequency  to  follow  a  triangular  distribution,  0  at 

x/d  =  0,  1  and  10^  sec  ^  at  x/d  =  .  5.  We  call  for  one  angle  of  incidence 

37.  5  deg,  and  one  frequency,  3  GHz.  We  ask  for  the  wave  to  be  incident 

from  the  right  and  polarized  in  the  TM  mode.  The  input  data  as  it  would 

be  entered  on  the  seven  cards  is  shown  in  Fig.  1.  The  corresponding 

output  is  shown  in  Fig.  2.  The  results  (in  this  case  the  real  and  imaginary 

parts  of  E  ,  B  ,  and  K,  and  the  value  ofVlE|^  +  primed  out  at 

y  X  X  ^ 

intervals  of  x/d  =  .  1  except  at  jumps  where  they  are  printed  out  twice, 
to  the  left  and  to  the  right  of  the  jump.  The  electric  field  is  always 
normalized  to  unity  on  the  incident  side.  The  fact  that  in  the  printout 
it  is  unity  at  x/d  -  1  serves  as  a  reminder  that  the  wave  was  required 
to  be  incident  on  this  side  of  the  slab.  On  the  transmitted  side  we 
have  a  single  plane  wave  in  free  space  for  which  E  /cB  =  -cos (37.  5°). 


-7- 


T 


9 

al 

I 


I 

S 


u 

o 

o 


u 


a  ^ 

a( 


I  I 


31 

!  * 

s 

s 

mm 

U 

9 


S 

8 

c 


if  ••!••••►•••►  ••♦Ml  ••  •* 

2  n  ««««»>  p«»  Mp>  •(>.  Am 

«i««««n«Aiiiw«n 

MMOOMOOMOMOM 

»  «  «n«>  «  oa^  aSS 

^O  O  *«P«K  MK  « 

2  •  A  lA #  o  om « 

2  ••••.•--*-*-m-«Vm*-*mV 

at 

M 

*  M  O  •  O  O  O  *«  p«  M  p«pp  ppppo 

^  M  MO  e  o  •  o  e  e e o  oe o 

M  «**«*«li|t|||* 

lylUUM  WMMMtuWktklkllil 

*  Op  5  ••••••►  o  ®  •«»»♦  o 

OKOO  OOOP  AMO  MOO 
g  •  OM  «  «  AmO  OOM  M»p  O 

*  •OO—'^AMO—OK  — AO 

2  A  AMO  «  0««M  A».M 

X 

^  ooo o ooo o ooo ooo 

O  «  "JO  o  ooo  o  ooo  ooo 

S,  u  WUM  Whllillil  WMM 

S  “I  o  P  M  A  O  MO  A  pppPM  OO  M 

*  ^  o  p p*n  p«p«o  M o o  o  oo  o 

1  S?;:S888S&&a::S$ 

,2  _  OK 

JO  •  » 

So  o  oooo222SS2**** 

•*-  I  ** 

2  422222r245P^^* 

.  TTTTT'J'?'?’?T'?T'?'f  *.•: 

223  ^ 

ooo  S  ooooooooSSooS  is 

•  ***  •  W  PPVWtfWW WWWW MW  I  ' 

•™-  I  gSISs^nSSSISS  -- 

-  8SSS^2:8SS8§8,  •• 

ooo  i  ----OOAOAAMMpSo  X 

o*o*’*p:  *  “ 

U  •^•«**»***p*p*p*«i*p***p«p«0  Mm 

K  ???T?????7??T?  St 

I  iiss^lHSH^i  r, 

m  A.;*'«v^oo*oVo*oV.:  *^2 

«  SS 

^  MX 

i  9  OOOOOOMAOOOO  3z 

S  H  •'M  A  O  A  O  O  OK  OO  O  0« 

^  oo  „  5^ 


•9- 


Figure  2.  An  Elxample  of  Output  Data 
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